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Buffet loads represent the primary source of high frequency loading for launch vehicles 
during the ascent portion of flight. Currently, experimental techniques establish the nature 
of buffeting using a rigid scale model of the vehicle. The buffet forcing functions resulting 
from such tests are then applied to reduced finite-element models of the full-scale vehicle to 
determine the response and consequent loading. This paper discusses the techniques 
required to translate model-derived, empirical buffet forcing functions into responses for the 
full-scale launch vehicle, as used to determine the buffet loading for NASA’s Ares I launch 
vehicle. 


I. Introduction 

B UFFET loads have been a key consideration for launch vehicles since the earliest years of the U.S. space 
program. Indeed, launch vehicles have been lost due to the high-frequency vibrations and loading that may 
occur due to buffet. Buffeting is closely associated with the transonic time of flight, where oscillating shocks and 
regions of unsteady, separated flow produce random aerodynamic inputs on a vehicle over a wide range of 
frequencies. This results in excitation of a variety of modes ranging from the vehicle’s primary bending modes to 
higher-frequency component and panel modes. 1 

In the present state of the art, buffet loads are predicted empirically via the measurement of scaled buffet loads 
on a model. The phenomena associated with buffet are too complicated to be accurately captured via CFD. Buffet 
is sensitive to the specific geometry of individual launch vehicles, with features such as geometric transitions and 
protuberances strongly influencing buffet loads. Therefore, in experimental testing it is critical that the geometry 
being tested accurately reflect that of the actual vehicle. 

“Hammerhead” configurations - such as that used by NASA’s Ares I vehicle - can lead to severe buffeting. 
Hammerhead vehicles feature a large-diameter forebody such as a payload shroud followed by a slender afterbody, 
usually the core of the launch vehicle itself Flow separation occurs where the body transitions from the enlarged 
front end to the smaller diameter, leading to buffeting. On Ares I, one such region is located close to the middle of 
the vehicle near the anti-node of the first bending mode. 2 These characteristics make buffet loads a significant 
consideration in the overall vehicle loads for Ares I. 

For the NASA Ares I launch vehicle, a rigid buffet test was conducted at NASA Langley's Transonic Dynamics 
Tunnel. 2 Here an instrumented, sting-mounted model of Ares I in its first-stage flight configuration was rigidly 
mounted in the wind tunnel. Flows over a range of velocities and dynamic pressures simulated the conditions 
experienced by the vehicle in the transonic and low supersonic flight regimes. 

In the 3.5% Ares I Rigid Buffet Model test, 256 pressure transducers monitored pressures at discrete locations 
over the surface of the model, with sensors concentrated in regions of geometric transition. These unsteady pressures 
were then integrated over effective areas to obtain the net unsteady force acting on the rigid buffet model in the axial 
and lateral directions. Via Strouhal scaling, these model-scale forcing functions were converted to full-scale 
equivalents for use in modeling the Ares I vehicle buffet loads. These scaling relationships can be found in a report 
on a rigid buffet test on a hammerhead version of the Atlas-Centaur I vehicle. 3 

The time-varying buffet forcing functions from NASA Langley were then applied to a flexibility model in order 
to determine the buffet response. For the Ares I buffet analysis, a model of the integrated vehicle was provided by 
NASA's Marshall Space Flight Center. This model, based on finite element representations, was reduced via a Craig- 
Bampton reduction. 
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To conduct a buffet analysis, buffet forcing functions are applied to the model’s degrees of freedom, and the 
forced response of the model is determined. In the case of this analysis, experimentally-derived buffet forcing 
functions were provided at discrete stations based on the geometry of wind tunnel model used; these stations 
differed from the force application points in the finite element model. This necessitated the development of a 
mapping between the two geometries. Several mapping methodologies were studied for the Ares I buffet analysis, 
and the resulting impacts of these mapping choices on the vehicle buffet loads are presented. 


II. Mapping Methodologies 


To determine the dynamic response of Ares I to 
buffeting, buffet forcing functions must be applied 
to the flexibility model's force application stations. 

These stations are physical locations 
(predominantly on the centerline of the vehicle) 
that have been retained in the model reduction p 
specifically as load-application points. In this case, r b 
the buffet forcing functions were provided at points 
taken from the location of sensors and pressure 
taps on the wind tunnel model. These points did 
not coincide with the force application stations of 
the flexibility model, so it was necessary to 
construct a mapping in order to apply the forcing 
functions to the model. 

In order to preserve the fidelity of the buffet 
forcing functions, it is important that the selected 
mapping conserve forces and moments. If the 
buffet loads were purely static, these criteria would 
be adequate. However, since the buffet forcing 
functions are provided in the form of time 
histories, it is also important to consider the phase 
relationships involved when two forcing functions 
are loaded at the same location. 

This section develops several different 
techniques for mapping buffet forcing functions to 
the finite element model. Section IV presents 
results for buffet loads obtained using these 
different methods. 



Figure 1. Two-point buffet forcing function mapping. A 

forcing function at point b can be mapped to flexibility model 
stations at a and c, provided that the force mapping 
conserves forces and moments. 


A. Two-Point Mapping 

The most basic mapping technique transfers loads to points on flexibility model and conserves forces and 
moments but makes does not preserve phasing in doing so. This method is exact for the static case, however. A 
buffet forcing function cannot be transferred to a single point on the Ares I model while still conserving moments, so 
it is necessary to apply the load to two points with appropriate weights so that the moment is balanced. 

Consider the schematic seen in Figure 1. A force, F b is applied at some generic Buffet Forcing Function 
(abbreviated here as BFF) station b. This force can be converted to a pair of forces applied at stations a and c that 
when summed equal the original force. The net moment applied by this pair of forces must be the same as that 
generated by F b . This will be the case if the net torque generated by F a and F c about point b is zero. 
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The following choices for F a and F c satisfy the force and moment conservation requirements. 




F , = Fdrrr) 

Xi+x 2 


This weighting satisfies the force summation requirement since 


F a +F c = F b (-^-) + F b (-^—) 
a c dk x 1 +x 2 *i+*2 


P a + Fc = h&^r) 


x 1 +x 2 ' 


Therefore, the two component forces sum to equal the magnitude and direction of the original: 


( 2 ) 

(3) 

(4) 

(5) 


Pa + Pc = P b (6) 

By choosing weighting factors based on the moment arms of the component forces, a net torque of zero is 
produced about point b. Taking counter-clockwise moments to be positive, the moment created by F a about point b 
is M a = —F a x 1? and that due to F c is M c = F c x 2 . The sum of the moments is therefore M b = F c x 2 — F a x 
Expanding the expressions for the point forces, it is shown that there is no net moment generated about point b. 

M b = P b (~~) x 2 ~ = 0 ( 7 ) 


Furthermore, it can be shown via a similar argument that the moment created by the mapped forces about an 
arbitrary point d is also conserved. The moment about d created by F b is M b = —F b (x). The combined moment due 
to F a and F c is M ac = —F a (x + x x ) — F c (x + x 2 ). Substituting the expressions for F a and F c , the moment can be 
written as 

Mac = -Pbi^Kx + *i) - P b (^- X ;)(* - * 2 ) (8) 

Combining the like terms yields 

(*1 + *2)(rg) = x 2 (x + x a ) + X r (x- x 2 ) (9) 

Expanding the distributed terms produces 

(*! + x 2 )(— ) = xx 2 + x x x 2 + xx x — x 1 x 2 (10) 

~ F b 

(*l+*2)(T^) = *(*!+ *2) (11) 

(^) = ^ (12) 

Finally, this result is identical to the expression for the moment about point d due to a single force at b. 
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B. Two-Point Mapping with Correlation Minimization 

While the two-point mapping ensures the conservation of forces and moments between the two sets of 
coordinates, the method makes no corrections for phase cancellations which may occur when multiple time-varying 
loads are combined at a single point. It was suggested that if multiple forcing functions were combined in a way 
that minimized the correlation between them, this would reduce the potential phase-cancelmg effects of this 
method. 4 

In general, phasing of the forcing functions presents the greatest concern when coupling occurs between the 
fluid and structural systems. However, in the case of the Ares I Rigid Buffet Model, it is expected that any structural 
interaction will be eliminated via the use of a rigid model. As a result, the correlation of the forcing functions along 
the length of the vehicle is already quite low, with correlation coefficients in the range of 0.1 or lower. 

The full-scale Ares I buffet forcing functions used in this analysis are each about 433 seconds in length. 

Previous temporal convergence studies have shown that only a 30-60 second window of each forcing function is 
needed to obtain sufficient convergence of the buffet response, and a 30-second transient provides very good 
convergence with substantially reduced computational overhead. 

As a result, it is possible to pick the ‘"windows” appropriately to minimize the correlation between summed 
forcing functions. AMATLAB program was designed that partitioned each buffet forcing function into 30-second 
windows, then picked forcing function blocks in a way that minimized the correlation. The first forcing function 
assigned any point is unique and can be chosen from an arbitrary window. When summing additional forcing 
functions, the program computes the correlation between the existing forcing function and all available 30-second 
windows for the new forcing function. The program then chooses the block that minimizes correlation. Using this 
method, the correlation coefficients were reduced by an order of magnitude to near zero. 


C. Four-Point Mapping 
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The preceding two-point mapping 
methodology transfers loads to pairs of 
points in the new coordinate frame. It was 
suggested that mapping the buffet forcing 
functions to greater numbers of points would 
yield more accurate results by expanding the 
range of influence of each forcing function, 
distributing the function over more stations 
on the finite element model. 4 A four-pomt 
mapping method was developed m order to 
validate this suggestion. 

In this method, each buffet forcing 
function in the source coordinate frame is 
mapped to four points in the new frame, 
bracketed by two points on each side as 
shown m Figure 2. 

The four-pomt method essentially splits 
force Fb into two forces, then uses the two- 
point method to apply part of the force at b 
to an inner pair. The remainder of the force 
is applied to a second, outer pair of points, 
again using the two-point method. Since the 
force and moment are conserved for each 
instance of the two-point method, then the 
total force and moment are also conserved in 
this implementation of the four-point 
method. 

This method is effectively a cascade of two-point mapping elements. In principle, more stages could be added to 
the cascade to further distribute buffet forcing functions. However, as the cascade becomes more complicated, it 
becomes increasingly difficult to trace the influence of an experimental forcing function on the stations of the finite 
element model. 


Figure 2. Four-point buffet forcing function mapping. A 

forcing function at point b can be mapped to four flexibility model 
stations instead of two. 


4 

American Institute of Aeronautics and Astronautics 


Although the four-point method has greater complexity than the two-point method, it is possible to use an 
identical correlation-minimizing algorithm. One possible advantage to this method is that each buffet forcing 
function can be assigned a larger region of influence than with the two-point method, and the weights of the “inner” 
and “outer” sets of forces can be tailored to achieve the desired force distribution. Note also that typically the 
weights applied to the mapped forces sum to unity; this ensures that forces and moments are conserved. If for some 
reason a non-conservative force distribution is desired, it is possible to use weights with a sum greater than one. The 
impacts of such weighting choices were studied, and these results are presented in Section IV of this paper. 


III. Buffet Analysis Methodology 

The Ares I buffet analysis made use of multiple finite element-based flexibility models corresponding to discrete 
times of flight. These models were similar structurally but differed principally in their mass distribution due to fuel 
bum. The entire phase of 1 st stage ascent was divided into a series of analysis periods, each utilizing a 
corresponding, stationary mass model. Previous studies have indicated the validity of this technique, provided that 
the periods are sufficiently short in duration. 5 

For a comprehensive buffet loads analysis, buffet would be modeled for a selection of representative trajectories 
and dynamic pressure (q) values. This technique is used to estimate the variation of flight loads with trajectory and 
flight conditions. However, this analysis of forcing function mappings studied only a single reference trajectory 
family. This approach allowed the impacts of the forcing function mapping methods to be assessed while keeping 
trajectory-dependent variables constant. Four analysis periods were created, spanning the times of flight where 
buffet loading typically occurs. The parameters assigned to each analysis period are shown in Table 1. 


Period 

Center Mach No. 

Condition 

Finite Element Model 

Buffet Forcing Function 

1 

0.7 

Transonic Flight 

t = 40 model 

M= 0.82, a= 4 deg 

2 

1.05 

Transonic Flight 

t= 40 model 

M= 0.98, a = 4 deg 

3 

1.4 

Max q-a 

t= 50 model 

M= 1.20, a = 4 deg 

4 

1.8 

Max q 

t=60 model 

M=1.55, a=4 deg 


Table 1: Description of run conditions for buffet mapping analysis. 

To determine the buffet loads for each analysis period, a 30-second window from the buffet forcing function was 
applied to the model. A modal transformation was applied to the model, and modes above 60 Hz were truncated. 
This truncation is a standard practice, because buffet forcing functions contain content to only about 50 Hz, and the 
higher vehicle modes in the finite element model cannot generally be trusted. High frequency buffet tends to excite 
components rather than the vehicle modes, and different modeling techniques are required to adequately capture 
buffet loading. 

The modal system was then integrated using a fourth order Runge-Kutta method. Following a transformation of 
the integrated response back to the original set of coordinates, the centerline forces and moments were recovered 
through the use of load transformation matrices (provided with the finite element models). The entire buffet analysis 
process was implemented as a series of MATLAB scripts. 

Note that when examining buffet loads, the key parameter is typically taken to be the RMS level of the unsteady 
buffet loading, scaled to the appropriate level of enclosure. Buffet loads follow a Rayleigh distribution rather than a 
normal (Gaussian) distribution; consequently, buffet loads are not reported in terms of the standard deviation. 

IV. Results 

The impacts of the different buffet forcing function mapping techniques may be observed through their effects 
on forcing function RMS levels at locations over the length of the Ares I vehicle. Figures 3-5 compare the centerline 
RMS force levels for the various mapping techniques with the nominal levels of the original, unmapped forcing 
functions as these exist in the wind tunnel model coordinates. A linear interpolation was performed to convert load 
levels in the wind tunnel model frame to the coordinates used by the finite element model. Only data from the 
M= 0.98, a = 4 degree case are shown, as trends for the other Mach numbers are very similar. 

Note the discrepancies observed between the interpolated original data set and the mapped forcing functions. An 
ideal mapping would preserve the magnitude of the original forces, resulting in a zero-percent difference between 
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the mapped force and the original. However, significant differences in load levels arise because the forcing functions 
are distributed to a much larger number of points in the new coordinate frame. Large differences are also seen when 
the original and mapped locations of forcing functions are far apart. These numerical issues, indicate that some 
discrepancy will inevitably exist between any mapped forcing function and its original counterpart. Therefore, the 
usefulness of Figures 3-5 lies in the relative trends and differences observed between mapping methods. 

For the two-point methods, little improvement is seen in results due to use of correlation reduction algorithms. 
Since the buffet forcing functions are relatively uncorrelated to begin with, this trend is not surprising. 

For the axial buffet forcing functions seen in Figure 3, it is clear that the 4-point uncorrelated method results in 
RMS levels typically lower than those of the 2-point methods. This is especially true for locations near the Ares I 
crew module, a region of major interest with regard to buffet loads. 

The impacts of the mapping methods on lateral forcing functions are shown in Figures 4 and 5. The distribution 
of lateral buffet forcing functions is relatively continuous along the vehicle, while the axial stations are concentrated 
only at regions of sharp geometry change, which contribute significantly to axial buffeting. As in the axial case, the 
4-point method is seen to result in lower loads than the 2-point methods. Use of the correlation reduction techniques 
does not greatly reduce the applied loads. Figure 4 shows the Y-direction results: the Z-direction results seen in 
Figure 5 have similar trends. 

Based on these results, it is observed that the two-point buffet mapping method provides a high-quality mapping 
as well as a straightforward implementation. It should be noted that these results are specific to buffet forcing 
functions which are largely uncorrelated. For highly correlated signals, special considerations must be made to 
preserve the phasing of the transients. In the case of the relatively uncorrelated buffet forcing functions analyzed 
here, the correlation-reducing methods were found to offer little performance improvement and did not have a 
significant impact on the buffet forcing function levels. 
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Figure 3. Effect of buffet mappings on axial RMS levels. 
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4. Effect of buffet mappings on lateral RMS levels, Y direction. 
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Figure 5. Effect of buffet mappings on lateral RMS levels, Z direction. 
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V. Summary 


This paper developed techniques for mapping buffet forcing functions from one model geometry to another. In 
the case of the Ares I buffet loads analysis, these methods were used to transfer buffet forc ing functions from a wind 
tunnel model coordinate frame to a finite element model. Analogous techniques could also be used to transfer loads 
from one finite element to another as models are periodically updated and improved. 

Two mapping techniques were developed: a simple, two-point mapping that preserved static force and moment 
relationships, and a four-point mapping also conserved these quantities while allowing greater tailoring of the buffet 
force distribution. These techniques do not conserve dynamic forces and moments because of the potential for phase 
cancellation when dynamic signals are added. Therefore, a technique was developed to minimize the correlation of 
additive signals, decreasing the likelihood of phase cancellation. 

In practice, the buffet forcing functions used in this analysis were found to be sufficiently uncorrelated to begin 
with; consequently, the correlation reduction method had little effect on buffet force magnitudes. The two-point 
buffet mapping method was seen to better preserve the character of the original buffet forces. The two-point method 
maps each buffet force to the minimum number of points needed in the new coordinate frame. 

The four-point mapping method led to a somewhat pathological effect: the multiple force distribution locations 
used in this method resulted in forcing functions dissimilar to their original counterparts. Forces were distributed 
further from their locations m the original coordinate frame. This issue could be corrected by using a non- 
conservative force weighting, however this technique was seen to increase buffet loads in a way that could not easily 
be quantified. 

This analysis focused on the effects of mapping on the applied buffet force RMS levels. Future work should also 
consider the impacts on the resulting buffet loads and accelerations, which are important considerations for any 
buffet analysis. In addition, mapping effects should be studied m the frequency domain, using the Power Spectral 
Density of forces, loads, and accelerations to assess any frequency impacts created by the mapping techniques. 
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